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ABSTRACT 

We consider the curvature radiation of the point-like charge moving relativistically along 
curved magnetic field lines through a pulsar magnetospheric electron-positron plasma. We 
demonstrate that the radiation power is largely suppressed as compared with the vacuum case, 
but still at a considerable level, high enough to explain the observed pulsar luminosities. The 
emitted radiation is polarized perpendicularly to the plane of the curved magnetic filed lines 
coincides with that of extraordinary waves, which can freely escape from the magnetospheric 
plasma. Our results strongly support the coherent curvature radiation by the spark-associated 
solitons as a plausible mechanism of pulsar radio emission. 

Subject headings: magnetic fields — pulsars: general — radiation mechanisms: nonthermal 

1. Introduction 

Although almost 35 years have passed since the discovery of pulsars, the mechanism of their radio 
emission still remains unknown. This is one of the most difficult problems of modern astrophysics. Soon 
after discovery, curvature radiation was suggested as a plausible mechanism for the observed pulsar radio 
emission (Radhakrishnan & Cooke 1969; Komesaroff 1970; Ruderman & Sutherland 1975). In fact, curvature 
radiation is the most natural and practically unavoidable emission process in pulsar magnetosphere. Most 
of the pulsar models suggest creation of dense electron-positron plasma near the polar cap. These charged 
particles move relativistically with Lorentz factors 7 ~ 10 2 — 10 3 along dipolar magnetic field lines. Therefore, 
they emit curvature radiation at the characteristic frequencies v c ~ 7 3 c/r c , where r c is the radius of curvature 
of field lines, which falls into the observed pulsar radio band if r c ~ 10 8 — 10 9 cm. This is a typical value of 
the radius of curvature of dipolar field lines at altitudes of about 10 7 — 10 8 cm, where the observed pulsar 
radio emission is supposed to originate (Kijak & Gil 1998, and references therein). 

However, some serious problems are encountered while considering the curvature radiation in pulsars. 
The first, and probably the most important one, is related with coherency of the pulsar radiation. It is 
well known that the incoherent sum of a single particle curvature radiation is not enough to explain a very 
high brightness temperature of pulsar radio emission. Therefore, one is forced to postulate the existence of 
charged bunches containing at least 10 15 electron charges in a small volume that can radiate the coherent 
curvature emission. However, it is not easy to form such charge bunches (see Melrose 1992, for review). 
Moreover, even if a bunch can be formed, it is not automatic that it will emit coherent radiation. For 
example, bunches formed naturally by linear electrostatic waves (e.g. Ruderman & Sutherland 1975) cannot 
provide any emission (Melikidze, Gil & Pataraya 2000, hereafter MGP00, and references therein). The 
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natural mechanism for the formation of charged bunches was first proposed by Karpman et al. (1975), who 
argued that the modulational instability in the turbulent plasma generates charged solitons, provided that 
species of different charge have different masses. Such charged solitons were observed in the laboratory 
electron-ion plasma (Sagdeev 1976) and perhaps even in the Earth ionosphere (Petviashvili 1976). In pulsar 
magnetospheric plasma, distribution functions of electrons and positrons are different, because plasma screens 
the electric field induced by co-rotation (Scharlemann 1974; Cheng & Ruderman 1977). This causes the 
effective relativistic masses of electrons and positrons to be different, which can result in a net charge of 
solitons formed in pulsar magnetosphere (Melikidze & Pataraya 1980, 1984, MGP00). The net soliton charge 
can also be induced by admixture of ions in the plasma flow above the polar cap (Cheng & Ruderman 1980; 
Gil, Melikidze & Geppert 2003). One should mention here that to explain coherent radio emission we do not 
necessarily need stable solitons but only large scale (as compared with Langmuir wavelength) charge density 
fluctuations. 

The second problem is related to the fact that the bunches (solitons) are surrounded by the magnetized 
plasma, which strongly affects the radiation process. The soliton size, which is determined by the level of 
the turbulence, is evidently larger than the wavelength of the Langmuir waves. A bunch is unable to emit 
radiation with the wavelength shorter than its longitudinal size. Therefore, the soliton should emit only at 
frequencies uj below the frequency of the plasma waves. In the pulsar frame of reference, the corresponding 
condition writes uj < 2y/jujp, where 7 is the Lorentz factor of plasma motion, uj p — (AiTe 2 n/m) 1 / 2 is the 
plasma frequency, n is the total number density of electrons and/or positrons, e is and m is the charge and 
the mass of electron, respectively. Since at the expected emission altitudes v p — uj p /2it ~ 1 GHz (see eq. [3] 
in MGP00) and 7 <~ 10 2 , the coherent radiation should be emitted at frequencies v< 20 GHz, as observed in 
radio pulsars. 

There are three waves propagating in the pulsar plasma: the extraordinary wave {uj ~ kc) polarized 
perpendicularly to the plane set by the ambient magnetic field and the wave vector k, and the two ordinary 
waves polarized in this plane, that is, the supcrluminal wave (uj > kc) and subluminal wave (uj < kc). Under 
the condition w < 2y/juj p , emission of the superluminal wave is heavily suppressed by the Razin effect (Razin 
1960). The subluminal wave may, in principle, be emitted, but it cannot escape from the plasma. In fact, it 
is ducted along the curved magnetic field lines preserving direction of the wave vector and eventually decays 
as a result of the Landau damping (Barnard & Arons 1986). The extraordinary wave escapes freely but 
in the infinitely strong, straight magnetic field it cannot be emitted because it does not interact with the 
plasma particles. However, in the curved magnetic field such wave can be emitted like in the vacuum case, 
in which a significant fraction of the curvature emission is polarized perpendicularly to the magnetic field 
line plane (e.g. Jackson 1975; Gil & Snakowski 1990). 

The aim of this paper is to calculate the curvature radiation within the plasma under condition 

£=-^«l. (1) 

We consider emission from a point charge (modelling a small soliton/bunch) moving circularly in the cylin- 
drically symmetric background. The magnetic field, which is assumed to be infinitely strong, has circular 
field lines, and the electron-positron plasma moves relativistically (with Lorentz factor 7) along these field 
lines. We find that the point charge emits curvature radiation in the extraordinary mode. The power of this 
radiation is suppressed by a factor s (eq. [1]), as compared with the vacuum case but still remains rather 
high. It was found recently that the Vela pulsar emits radio waves polarized predominantly in the direction 
perpendicular to the plane of dipolar magnetic field lines (Lai ct al. 2001), strongly suggesting that these 
are escaping extraordinary waves. Our results demonstrate that such radiation can be generated within 
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the pulsar magneto-sphere by means of the coherent curvature radiation. This radiation can also leave the 
magnetospheric plasma and reach a distant observer. 

MPGOO proposed the spark-associated soliton model for pulsar radio emission, in which the charged 
solitons are formed by nonlinear evolution of plasma waves triggered by a sparking discharge of a high 
accelerating potential drop above the polar cap. These relativistic solitons emit the coherent curvature 
radio emission at altitudes of about 50 stellar radii in typical pulsars. The total luminosities emitted by an 
ensamble of solitons are consistent with the observed pulsar fluxes. The formation of charged solitons in 
pulsar magnctosphere, as well as their coherent curvature radiation, is well justified by MGP00. The only 
deficiency of this model is that the curvature radiation was calculated in vacuum approximation. We argue 
in this paper that the soliton coherent curvature radiation is still a plausible mechanism of the observed 
pulsar radio emission even when the presence of the ambient plasma is taken into account. 



2. Wave equations 

We now examine the curvature radiation from the point-like charge moving relativistically along a curved 
trajectory within the electron-positron plasma. We find the radiation power by calculating the work done 
on this charge by the excited electromagnetic field (e.g. Melrose 1980). It is convenient to considered the 
curvature radiation in the cylindrical coordinates (r, 9, z). The force lines of the infinitely strong magnetic 
field are assumed to be concentric arcs with the radius of curvature r being much larger than the size of 
the considered region of emission, marked by dotted parallel lines in Figure 1. The charged particle moves 
along the dashed circle with radius R. The thick shadowed line represents the actual magnetic field lines, 
which are locally modelled by concentric arcs. The thick dot represents the pulsar and r a is a radio emission 
altitude. 

The wave equation for the electromagnetic vector potentials writes 

1 d 2 A . . 4tt 



^ A = ~ J ' (2) 
provided that vector and scalar potentials are subject to the Lorentz gauge condition 

1 a* 

-- + V.A = 0. (3) 

The scalar potential ^ satisfies the wave equation as well, provided that the total charge of the system is 
conserved. The electric field of the wave is defined as 

■ ~i£-vt. (4) 

We use equations (2) and (3) as a complete system of equations. Let us note that in the infinitely strong 
magnetic field the particles are forced to move exclusively along the field lines and, consequently, the current 
can be excited only along the ^-direction. 

Making the Fourier transform 

f(r,s,k z ,uj) — J exp(iujt — is8 — ik z z)f(r,6, z,t)dtd6dz, (5) 

where / denotes it and components of A, J and E, we get the following set of equations 

* = a ±A e + ^A z -i±£-rA r ; (6) 
ur uo lot Or 
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1 d A fto 2 s 2 ,,\ , „ s A 
^ 1 d A fu 2 s 2 , 2 \ „ „ s „ 4tt t 

— -— rA e + — - -j - k 2 z )Ag + 2i—A r = J ; 8 

or r or \ & r z j r z c 

AlA r A (^L- S 1-]A A i) (9) 
dr r dr z \c 2 r 2 z ) z 

One should remember that according to the Landau prescription, ui is assumed to have a small positive 
imaginary part. We do not write down it explicitly here but will restore it in the appropriate place (sec 
bellow eq. [24]). 

Now we need to find the current density associated with the point charge JP art and current density 
jpiasma excited in plasma by the wave. If we assume that the emitting particle rotates with the velocity V 
(the corresponding Lorentz factor is T) along a circle with the radius R, we can introduce the corresponding 
current density as jg — (QV/R)5(r — R)S(0 — Vt/R)S(z), where Q is a charge of the emitting particle. Making 
the Fourier transform we get 

J part =2nQ^6(r-R)6(w-s^y (10) 

In the infinitely strong magnetic field, only 0-component of the electric field of the emitted wave excites the 
currents in the plasma. Assuming that the plasma is cold and moves with the velocity v, we get 

jpla S ma = .^l W 

4^ 7 3 (w _ s h) 2 

(see details in Appendix A). Here 7 = (1 — v 2 jc 2 )^ 1 ! 2 in the plasma Lorentz factor. The Fourier component 
of the longitudinal electric field (eq. [4]) can be expressed via the potentials in the form 

Eg=i{^-A 6 -Ul). (12) 

Our objective is to find such solution of equations (6-12), which remains finite as r — > and comprises an 
outgoing wave as r — > 00. 



3. Short wavelength approximation 

Ultrarelativistic particles emit waves with the wavelength much shorter than the curvature radius of 
their trajectory. Therefore, we can solve the wave equations in the short wave approximation, uir/c ~ s 1. 
Let us start with the so-called WKB approximation, in which A cx exp(i J k r dr) and therefore the solution 
splits into two waves 

u? s 2 
c 2 

and 



k r - ~J 12 ~ k z ( 13 ) 



- 2 -(S + *= + *O c "w^if(" 2 -^ 2 )' (14) 

Equation (13) describes the extraordinary wave (u> — kc) polarized perpendicularly to the plane set by the 
local magnetic field and k vector, while the ordinary wave polarized in this plane is described by equation 
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(14). Under the condition expressed by equation (1), equation (14) has a solution corresponding to the 
nearly transverse wave with the dispersion law in the form 

- 2 - (^) 2 (i - (fc y| c2 ) if + ® ° 2 « < (i5) 

The WKB approximation is not valid close to the classical reflection radius r , defined as 



U? - - k'tc Z = 0. (!(,) 



At r w ro, the refraction indices AT = fcc/w of extraordinary (eq. [13]) and ordinary (eq. [15]) waves are very 
close to each other (uj — kc), thus providing conditions for the wave coupling. In this region k r goes to zero, 
thus forming a caustic zone. 

Consequently, we should solve differential equations (6-12) for the electromagnetic potentials in the 
region r w r and then match the obtained solution with the appropriate WKB solution (Eqs. [13] or [14]). 
Note that because the particle emits nearly along the local magnetic field, the radius of the particle trajectory 
R should be also close to r . It is therefore convenient to use the radial variable 

- r - r V/3. (17) 



In the region of interest x <~ 1 and correspondingly \r — ro\ ~ ros~ 2 / 3 <C ro. The power series expansion in 
small s -1 / 3 <C 1 allows to reduce equations (7-9) to the form 

A" + 2xA r - 2is- 1/3 A e = 0, (18) 

Ag + 2xA e - 2is- 1/3 A r = D, (19) 
A" + 2xA z = 0, (20) 

(see Appendix B for details), where 

D = -47rr2 s - 4 / 3 (jP lasma + jP art ) 

y p7 (2(x + a)A e + is^A' r ) q „ 2t ^„_ 2/3 ,^ (2 1) 



(l + 2 S - 2 / 3 7 2 (a + x)) 



,■,1/3 4' / 

r} 2 - 8n 2 VQs- 2 / 3 S (x — X)5\u 
a + x)) V 



R 



1 s 2 ' 3 , 1 s 2 / 3 k 2 c 2 



X = -a--— and a = -— -f-, (22) 

where prime denotes differentiation with respect to x-coordinate and X denotes x-coordinate of the emitting 
particle (see Appendix B) . We do not assume any relationship between A r and A g and thus we retain terms 
proportional to a small factor s -1 / 3 in equations (18) and (19). 

The gauge condition expressed by equation (3) still permits a gauge transformation $ — > \& — d$/dt, 
A — > A + V$, provided that $ satisfies d'Alembert's equation, which reduces to 4>" + 2x$ = if lot ^> 1 
and |r — r | <C r - Therefore, one can choose A z — (c.f. eq. [20]) and solve equations (18) and (19) together 
with equation (3) or equation (6) . The solution must comprise an evanescent wave at x < X and an outgoing 
wave at x > X. 
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The Green function for equations (18) and (19) can be found from the solution of the corresponding 
homogeneous equations matched at the point x = X according to the following matching conditions 

A' g (X + 0) - A' g (X - 0) = -8n 2 QVs- 2 / 3 S (u - s-jQ and A' r (X + 0) = A' r (X - 0), (23) 

which can be obtained straightforwardly integrating equations (18) and (19) between limits X — and X + 0. 

Let us consider the first term of D (see eq. [21], which is just the right hand side of eq. [19]). The 
nominator in this term (which according to eq. [B9] is proportional to Eg) should be small, due to the 
condition expressed by equation (1). Thus, in the zeroth approximation in e (eq. [1]) the following relationship 
holds 

1/3 

2(x + a + i0) r x ' 

Here we take into account the fact that according to the Landau prescription u and consequently a have a 
small imaginary part. Substituting the obtained relation into equation (18) we get a closed equation for A r . 
Let us note, that equation (24) means that Eg — 0. Therefore in order to calculate the emission power we 
need the next approximation in e (see the next section). 

Now it is convenient to introduce the dimcnsionless function u{x) such that 

8n 2 QV X + a 

A r = t^iSOjj - sV/R)u, (25) 

s X 

(see Appendix B for details), which satisfies the following equation 

u' 

u" + 2xu = 0. (26) 

x + a + iO 

The matching conditions (eq. [23]) can be now reduced to 

u(X + 0) - u(X - 0) = 1 and u'(X + 0) - u'{X - 0) = 0. (27) 

Equation (26) has the Airy type asymptotics u — ► exp(±(2v / 2/3)«a; 3 / 2 ). Therefore, we should look for the 
solution which satisfies the conditions u — > cxp(— (2\/2/3)(— x) 3 / 2 ) asu — oo and u — > exp(i(2-\/2/3)x 3 / 2 ) 
as x — > oo. Then we can find A r and Ag from equations (24) and (25), and match them to the WKB solution 
(eq. [13]) as x — > oo. 

It is straightforward to calculate the electric field of the waves in the far zone using equations (4), (6) 
and (24) as 

(E r ,E e ,E z ) oc 0, -^T) exp ^^x 3 A . (28) 



This solution is valid in the region 1 <C x <C s 2 ^ 3 (see eq. [17]), where it is smoothly matched to the global 
WKB solution expressed by equation (13). One can see that the outgoing wave constitutes the extraordinary 
mode polarized predominantly in the direction perpendicular to the plane of the curved magnetic field lines. 

Unfortunately equation (26) cannot be solved by using any known special functions, and has to be 
solved numerically. For the numerical solution, it is convenient to reduce equation (26) to the Rlkkati type 
equation 

w' + w 2 + 2x = 0, (29) 

x + a + iO 
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where w — v! /u. We need only the imaginary part of u, because, as it is shown in the next section, the 
emissivity of the particle is expressed via lmu(X) (see eq. [31,32]). Using the matching conditions expressed 
by equation (27), we can express lmu(X) by the solution of equation (29), that is 

lmu{X) = ^1^+ (30) 

v ' (w- -Rcw + ) 2 + (Imw+) 2 V ; 

Starting from a large negative x for which w = \/2\/—x (an evanescent wave) and integrating numerically to 
x = X, we found a solution w_ = w(X — iO). On the other hand starting from a large positive x, for which 
w = i\f2^px (an outgoing wave) and integrating to some x = e <C 1, then integrating (in the complex plane) 
along the semicircle x = ecxpicf), and finally integrating further along the real axis to the point x = X, wc 
found another solution w + = w(X + iO). 



4. Emissivity of curvature radiation 

The emitted power can be found by calculating the work done by the radiative electric field on the 
emitting particle 

P = -J r)£7 fl (t,r) (31) 

(e.g. Melrose 1980). Note that Eq = in the zeroth approximation in e (sec eq. [24]), therefore we should 
find the next approximation. It can be done easily substituting into the left hand side of equation (19) the 
values found in the zeroth approximation. We find that 

Ee = 2^V 1 r 4 %g / 1 ' ^ 5 (s 4) u (*) • ( 32 ) 

Now we can write the emissivity in the form 

- ^ =£^-2!)%,,,, (33, 



' dtdvdfi it R 2 V T 2 

where 



2 2 

F ^^= (TT^) Imu(x) ' (34) 

and dimensionless frequency v and angle /x are defined as 

w R k z k z c 

v = — = cj — — and u = —1 = 1 . (35) 

lu c cT A ' k LU 

Let us note that a position X of the emitting particle can be expressed by means of v and fi as X = 
-0.5^ 2/3 (1 + ii 2 ). The parameter 



1 / cT 3 



2 



A = — — - < 1 (36) 

describes the suppression of the emissivity of the curvature radiation by the surrounding plasma. Actually, 
this is the ratio of the curvature radiation frequency to the plasma waves frequency, thus it corresponds to 
the condition expressed by equation (1) at the frequency ui = ui c = cT 3 / R (see also discussion below eq. [42]). 
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As a consistency check let us note that it follows from equation (33) that 77 = if V = 7 (no radiation from 
uniform plasma motion), as expected. Function F {v, \i) is plotted in Figures 2 and 3 for various values of v 
and \x. Let us note that \i = 1 corresponds to the radiation angle 1/T from the plane of the source motion and 
v = 1 corresponds to the characteristic frequency of curvature radiation v c = cT 3 /2ttR. The polarization 
of this radiation is that of the extraordinary wave. The numerical values on the vertical axes in Figs. 2 
and 3 can be useful in estimating the emissivity of curvature radiation of the point-like charge Q moving 
relativistically (Lorentz factor T) through the relativistic plasma (Lorentz factor 7) in pulsar magnetosphere. 
Note that the emissivity is zero if /x = 0, which means that although all the radiation is concentrated within 
the angle <~ 1/r around the plane of the charge trajectory, there is no radiation exactly in this plane. 

The total luminosity from the charge can be obtained as Li = J rjdfidi>, where rj is the emissivity 
described by equation (33). One can find by numerical integration that J F(/j, v)d\idv = 0.53 and obtain 



According to our results the escaping waves are polarized almost perpendicularly to the plane of the 
magnetic field curvature, whereas the wave itself propagates nearly in this plane. Thus, this wave constitutes 
the extraordinary plasma mode. At first sight this conclusion seems counterintuitive. It is generally believed 
that in the infinitely strong magnetic field the extraordinary wave cannot be excited because the electric 
field of this wave is perpendicular to the external magnetic field and therefore it does not interact with the 
plasma. However, this statement is correct only in the case of the straight magnetic field. It is the curvature 
of field lines which makes a difference. For example (Luo 1993) and (Luo & Melrose 1995) demonstrated 
that the extraordinary mode can be generated in the infinitely strong magnetic field provided that the field 
lines are twisted. In this paper we show that the extraordinary mode can be excited even if the curved field 
lines lie in the plane, as it is in the case of a dipolar pulsar magnetic field. 

To illustrate that let us consider the vacuum case first. The curvature radiation is polarized predom- 
inantly in the plane of the particle trajectory, however about 15% of the total power is emitted with the 
polarization perpendicular to this plane (Jackson 1975, p. 675). Of course, the radiation which is emitted 
exactly in the plane of the particle trajectory is completely polarized in this plane. However, the radiation 
emitted at the angle about 1/r to the plane of the particle trajectory (k z /k ~ 1/r) is elliptically polarized. 

In the case of the plasma in the infinitely strong magnetic field, a charged particle can emit only 
waves with a nonzero component of the electric field along the external magnetic field. According to the 
standard classification, these waves are referred to as the ordinary plasma modes. Strictly speaking the WKB 
approximation is violated in the radiation formation region, and therefore one cannot introduce the normal 
modes there. However, one can speak in terms of the normal modes heuristically. At the condition (1) the 
particle emits subluminal, quasi-transverse Alfven waves (the emission of the superluminal ordinary mode is 
suppressed by the Razin effect). The Alfven wave polarized in the plane of the curved external magnetic field 
does not propagate outwards, because it is ducted along the curved magnetic field and eventually decays 
by the Landau damping (Barnard & Arons 1986). However, just like in the vacuum case, there are waves 
with k z /k ~ 1/r and they are polarized perpendicularly to the plane of the external magnetic field line. 
The polarization plane of such a wave would rotate by the angle about unity as it propagates the distance 




(37) 



5. How can the particle emit radiation polarized perpendicularly to the plane of its 

trajectory? 
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I ~ Rk z /k ~ i?/r, if the adiabatical walking condition 



(iV-l)fcZ»l 



(38) 



was satisfied (Cheng & Ruderman 1979). Here iV is the refraction index of the wave, which in our case can 
be estimated as 



(see e.g. Barnard & Arons 1986, or eq.(15)). Taking into account that k ~ T 3 /i?, one gets straightforwardly 
(N—l)kl <~ 1. So, the adiabatical walking condition (38) is not satisfied and therefore the wave escapes from 
the plasma retaining the initial polarization in the direction perpendicular to the magnetic field lines plane. 
In the far zone this wave should be classified as the extraordinary wave. Thus the extraordinary mode can 
be emitted in the curved infinitely strong magnetic field via the linear coupling of the normal modes in the 
radiation formation region. 

It is interesting to compare the expression for the luminosity of the curvature radiation in plasma 
(Eq.[37]) with the power of single-particle curvature radiation in vacuum p = (2/3)cQ 2 T 4 R~ 2 (Jackson 
1975). As one can see, equation (37) has a quite intuitive form. The coefficient 0.1 ss (2/3)(l/7) is related 
to the fact that this luminosity corresponds to the "perpendicularly" polarized radiation, which contains 
only about 1/7 of the full power of curvature radiation (Jackson 1975). Thus, the term A (l — 7 2 /T 2 ) in 
equation (37) represents the influence of the plasma, while the factor 0.1 represents the weaker extraordinary 
mode of curvature radiation. 



The two-stream instability induced by a sparking discharge of the high potential drop above the polar 
cap easily excites the strong Langmuir turbulence at altitudes 10 7 — 10 8 cm (Usov 1987; Assco & Melikidze 
1998), where the pulsar radiation is supposed to originate (e.g. Kijak & Gil 1998). As shown by MGP00 this 
instability can lead to the formation of charged rclativistic solitons in pulsar magnctosphcric plasma, which 
may play a role of point-like charged bunches that can emit coherent curvature radiation. The size of the 
solitons is about the characteristic instability length-scale, which can be estimated as 



where W is the energy density of the Langmuir waves and the characteristic plasma "temperature" (the 
kinetic energy spread) is about mc 2 . All primed values refer to the proper plasma frame of reference, in 
which the transverse size of the soliton should be about the same as the longitudinal size (by the causality 
argument). Taking into account that in the pulsar frame the longitudinal size decreases 7 times whereas the 
transverse size remains the same, one can assume that our bunches have a form of pancakes. The soliton 
emits at wavelengths about its longitudinal size A ~ d = d' /j, therefore its transverse size is about -yd. Note 
that emission from two points separated by the transverse distance d± is added coherently only within the 
angle <~ X/d±. Because the relativistic emission is concentrated within the angle ~ I/7, one can see that 
the transverse size of a coherently radiating bunch cannot exceed <~ -fd. 




(39) 



6. Implications for pulsars and discussion 




(40) 



The solitons in a real turbulence do not survive for a long time. The stronger the turbulence, the shorter 
the soliton life time. The soliton emits as a single entity if it survives for a radiation formation time, which 
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is defined as a time during which the phase of the emitted waves changes by about 2tt in the emitter frame 
of reference. This time is about X/(cT) ~ d/(cT) in the soliton frame. Because the soliton cannot decay for 
a time shorter than it takes for the light to propagate across the soliton itself, the above condition does not 
place any restrictions on the proposed emission mechanism. What we need in fact are not stable solitons, 
but only long wavelength fluctuations, whose existence is guaranteed by the modulation instability criterium 
(MPGOO, eqs. [A20] and [A23]). 

The coherency condition requires that the radiated waveleght is much larger than the longitudinal size 
of the emitting soliton. This is equivalent to the condition expressed by equation (1). Adopting lo c = cT 3 /R 

1/2 / . \ 1 /2 

and ojp — (47re 2 n/m) , where n = kugj = 1-4 x 10 n Krg~ (P_i 5 /P\ cm~ 3 is the plasma number density 

at the altitude r = r 6 x 10 6 cm, ncj is the Goldreich & Julian (1969) number density, R ~ lOr is the radius 
of curvature of dipolar field lines, T = r 2 100, 7 = 72IOO and k is the Sturrock (1971) multiplication factor, 
we obtain 




which should be much less than unity (say at least 0.1), where P is the pulsar period in seconds and 
P_i5 = P/10~ 15 is the period derivative in units of 10 -15 s/s. The actual value of k is somewhat uncertain, 
with estimates ranging from 1 to 10 . The lower values of k < 100 correspond to the stationary acceleration 
models (e.g. Hibschman & Arons 2001), while in the non-stationary sparking scenario being of interest here, 
one can expect larger values of k « 10 2 — 10 4 (Ruderman & Sutherland 1975; Melrose 2000). However we 
will estimate e as a function of k treated as a free parameter. One can see from the above equation that for 
emission altitudes 10 < r 6 < 100 (Kijak & Gil 1998) e< 0.1 if k> 100 (for P = P_i 5 = 1, 72 = 2 and T 2 = 4). 
Thus, the basic condition (eq.[l]) for the spark associated soliton curvature radiation is certainly satisfied. 
Now let us consider the energetics of relevant processes. 



6.1. Radiation luminosity of the turbulent plasma 

Let us consider the emitting region within a tube of open dipolar field lines, whose radial extent Ar <~ r a 
(see Fig. 1) and the cross-section surface area is S. The total number of solitons can be estimated as 
N <~ Sr a /(2V), where V is the soliton volume, and we assumed that A 3> cir/y^yWp, thus e ~ A « 1 (eqs. 
[1] and [41]). Assuming reasonably that r ~ 2j (see MPGOO) the total luminosity of Af solitons may be 
estimated as 

L w 0.1AcQ 2 r 4 i?- 2 AA, (42) 

where Q is the point charge moving relativistically (with Lorentz factor T) along a circle with radius R and 
A = w 2 /4w 2 7 is the suppression factor defined by equation (36). The charge of the soliton is Q = eSnV, 
where Sn is the difference between the densities of positrons and electrons within the soliton. Since the 
charge density within the soliton is determined by the Langmuir wave energy density W (e.g. MGP00), then 
we have 

^~PV) 2 . (43) 
n \nmr7 / 

Noting that Q 2 M = (l/2)e 2 {5n/nf n 2 r a VS = (1/8tt) (I/7C 3 ) (Sn/n) 2 u 2 p r a VL , where L = nmc 3 1 S is the 
total kinematic power of the plasma flow (section 6.2), and remembering that the soliton volume is estimated 
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from equation (40), the total soliton luminosity can be written in the form 

V72/ \irnnc 2 ) 

Here we assumed that r a /R <~ 0.1 for dipolar field lines at altitudes r a ~ (10 7 — 10 s ) cm 
can be compared with the observed pulsar radio luminosity (sections 6.3 and 6.4). 

6.2. Pulsar kinematic luminosity 

Since the kinematic flux is conserved along the tube of open dipolar field lines, we can calculate L = 
nmc 3 -fS near the polar cap, where r 6 is about 1 and S = 7rr 2 ~3x 10 8 P _1 cm 2 is the polar cap surface 
area. Thus, the kinematic luminosity 

L « 10 2 V^-i5^ 3/2 erg a" 1 . (45) 

Since the product 0.05< P^i 5 P~ 3 ^ 2 < 10 2 , with a typical value being about unity and 7 <~ 100 we have a 
typical value £q ~ k10 26 erg s _1 , where n > 100, 



6.3. Observed pulsar radio luminosity 

The pulsar radio luminosities Lr can be obtained from measured fluxes and estimated pulsar distances. 
If I = log(S4ooD 2 )i where Szioo is the mean flux density at 400 MHz given in mJy and D is the pulsar distance 
in kpc, then 

Lr = 3.5 x 10 25+i erg s" 1 . (46) 

According to Figure 9 in Taylor, Manchester & Lyne (1993, Pulsar Catalogue), 1< l< 3, with median value 
< I >^ 2. Thus, the pulsar radio luminosities 3.5 x 10 26 < Lr< 3.5 x 10 28 (erg s _1 ), with median value 
C R ~ 3.5 x 10 27 erg s" 1 . 

6.4. Comparison with observed luminosities 

Now, using equations (41) and (45) we can rewrite equation (44) in the form 




Assuming 72 ~ 2, T 2 ~ 4 and k ~ 100 we can estimate the luminosity as L <~ 5 x 10 27 erg s _1 (close to the 
median value jCr), provided that (W/jm e nc 2 ^ ~ 0.5. This requirement does not seem too excessive. In fact, 
the two-stream instability due to overlapping of adjacent plasma clouds associated with successive sparks 
(Usov 1987; Asseo & Melikidze 1998) should provide a high enough level of turbulence, because energies of 
the plasma and the beam are of the same order. Thus, one can expect that about 50% of the beam energy 
will be transferred to the plasma waves. Let us note that the luminosity is quite sensitive to estimations of 
the soliton volume. The longitudinal size of the soliton is defined by the Langmuir turbulence (eq. [40] ) , but 
the cross-section may can be estimated using coherency conditions according to the radiated wavelength. 
Then the luminosity can be higher by at least one order of magnitude. Moreover, the term involving 72 and 



(44) 

This luminosity 
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T2 can drastically increase the value of L under minor changes of Lorentz factors. Therefore, the model of 
curvature radiation of the spark-associated solitons developed by MGPOO in vacuum approximation is still 
a plausible explanation of pulsar coherent radiation, once the influence of the ambient plasma is taken into 
account. The power of curvature radiation is suppressed by plasma but not drastically with respect to the 
vacuum case. The polarization of curvature radiation emitted in plasma is that of the extraordinary mode, 
so it can escape from the pulsar magnctosphere. 

Interestingly, Lai et al. (2001) argued recently that the polarization direction of radio waves received from 
the Vela pulsar is perpendicular to the planes of dipolar magnetic field lines. It is instructive to follow their 
argument, which strongly implies that the observed radiation from this pulsar represents the extraordinary 
plasma mode. In fact, Lai et al. (2001) were able to demonstrate that in the fiducial phase corresponding to 
the fiducial plane (containing the rotation and the magnetic axes as well as the line-of-sight) the radiation 
is polarized perpendicularly to the plane of the dipolar magnetic field lines. This argument can be extended 
to every phase within the pulse window, since the position angle swing in this pulsar is known to satisfy 
perfectly the geometrical rotation vector model (Radhakrishnan & Cooke 1969). This means that in the 
Vela pulsar the polarization of observed radio emission is consistent with the curvature radiation originating 
in pulsar magnetospheric plasma. 

This paper is supported in part by the Grant 2 P03D 008 19 of the Polish State Committee for Scientific 
Research. We thank Don Melrose for his hospitality during our staying at the School of Physics (Sydney 
University), where the final version of the paper was completed. We thank E. Gil and U. Maciejewska for 
technical help. 

A. Current in the plasma 

In order to calculate the current density in the plasma let us start with the equation of motion which 
in the relativistic case writes 

m 7 3 -^ = eEg. (Al) 

In the Fourier components we have 

dv ( d6\ 
— = — i \u — s— v. 
dt V dt J 

Since ddjdt — v/r, thus the equation of motion writes 



v=-^ ( Ee v (A2) 
7 d m (lu — ^v) 

The electric field is expressed by equation (12). Let us use the continuity equation in the form 

dn „ , . 

— + V • {nv) = 0. 

Assuming that the perturbation of the density can be expressed as n — > n + Sn, then for the corresponding 
Fourier components we can write dn/dt = —iutSn and V • (nv) = vWn + nV • v = iv(s/r)Sn + iv(s/r)n. 
Then using equation (A2) we obtain 

Sn=^ S — ^L_. (A3) 
7 3 m r (u-i v ) 
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Assuming that the unperturbed charge density equals to zero, we can express the charge density perturbation 
as p — 2eSn. Now we can find the current density from the following equation 

which in the Fourier transforms writes 

-iup + i-J e = 0. (A4) 
r 

Substituting equation (A3) into (A4) we arrive at equation (11). 



B. Wave equation in short wave approximation 

We are looking for the solution of equations (6 - 9) in an narrow region r — ro <C ro, therefore we can 
substitute r with ro into coefficients, which vary slightly in this region. Then we can expand the coefficients 
which vary significantly to the first order in r — ro, those are the coefficients in the curly brackets in equations 
(7 - 9), which become zero at r — r and the coefficient in the denominator of equation (11), which is close 
to zero at v — > c and k z <C co/c. Using the approximation s ^> 1, k z c <C w and 7 ^> 1 from definitions 
expressed by equations (16) and (17) we get 

s ( 1 r, 2 

U! = — C 



ro'V ' 2s 2 

and 



1 _ 1 

r r 

Then using a definition 



i + 5^J (Bi) 

- = -(l- s- 2 ^x] . (B2) 
r r n V / 



it is straightforward to obtain the following formulae 



, ,2 „2 4/3 

^ - % - k l = 2 ~ x (B4) 



and 




/ jk z c\ 2 _ 2/3 2 \ 1 sc 



, + 2s- 2 / 3 j 2 x = -L— (l + 2 S - 2 / 3 7 2 (a + xj), (B5) 



We also need an evaluation of expressions for the currents (see eqs. [10] and [11]). The argument of 
the first 5-function in equation (10) can be evaluated using the condition ujR = sV, which follows from the 
argument of the second ^-function. Therefore, we get 

2/3 

S(r-R) = - — S (x - X) , (B6) 

where 
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describes the position of the emitting particle. Using equation (B6), we can express the current density 
excited by the emitting particle (eq. [10]) as 



47rr-2 s - 4 / 3 J e part = S^s-^VQS {x-X)s(u>- s^j . 



(B8) 



Using equation (6) we can rewrite the expression for the electric field (eq. [12]) in the form 

(4/3 \ S. 

2 s x + k 2 z )A e -- S -^A' r , (B9) 

and with the use of equations (B5) and (B9) the expression for the plasma current density (eq. [11]) can be 
written as 

-47rrg s - 4 / 3 J 2 plasma = 8^ (l + 2s" 2 / 3 7 2 (a + xj) ^ (2 (x + a) A e + is^A' r ) . (BIO) 

Now using equations (B8) and (BIO) we obtain the right hand side of equation (19) as it is defined by 
equation (21). 

Now we need to find the matching condition for A r . Using equation (24) we have 

% s 1 / 3 1 
A' B = — : tK-i r^e, (Bll) 

and getting A" from equation (18) and substituting it into equation (Bll) we get 

A r = -i s -V 3 ^±^A' e . (B12) 

x 

Now using the matching condition for A' g (eq. [23]) we obtain the following matching condition for A r 

V\ X + a 



A r (X + 0) - A r (X-0) = i^VQs^S (u> - s^j 



X 



(B13) 



With the help of the above equation, the matching condition expressed by equation (27) follows from defi- 
nition expressed by equation (25). 



C. Emissivity 

The reverse Fourier transform writes (c.f. eq. [5]) 

00 00 00 

f (9,z,t) = --^y J J J f(s,k z ,u)exp(-i(u)t-k z z-s6))(Lidk z ds. (CI) 

—00 —00 —00 

Substituting the Fourier transforms of the current density and the electric field to equation (31) we obtain 



00 

P = g / rdrd9dzdtJg art (u/,r, s' ,k' z ) Eg (u>,r, s,k z ) 

(2n) J 



(27T) 

— uu 

x exp (i (u - u') t - i (k z - k' z ) z - i (s - s') 6) dwdk z dsdw' dk z ds' , (C2) 
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where Eg is a complex conjugate of Eg. Then using the following relation 

oo 

J exp {% (k z -k' z )z + i{s- a') 9} dOdz = (2nf S (k z - k' z ) S {s - s') , (C3) 

— OO 

we obtain the emissivity of the charged particle in the form 

oo 

d^k = J rdrJ ^ ^'^ E e («,r )fll *)e-^)W. (C4) 

— OO 

Now we can integrate over s and uj' using 

oo 

J s(u' -s^S^s-GJ^e-^-^dsduj' = 1, (C5) 

— oo 

and over r using delta-function 6(x — X) in the expression for Jg art (see eq. [10]). Then substituting uj and 
k z by the dimensionless values v and fi (eq. [35]) we obtain equation (33). 
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Fig. 1. — Schematic picture of the radiation region (marked by dotted parallel lines) and the system of 
cylindrical coordinates (r, 0, z), with z-axis directed upright to the plane of the figure. The charged particle 
moves along the dashed circle with radius R, which is locally tangent to a dipolar magnetic field line marked 
by a thick shadowed line. The thick black dot represents the pulsar and r a is radio emission altitude. 
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Fig. 2. — Distribution of the radiation power emitted in all directions F(v) = j F(v, /j,)d[i (lower panel) and 
in various directions F(v,m) (upper panel) over the dimensionless frequency v = u)/ui c = ljRc~ 1 T~ 3 , where 
/ii = 0.5, /Lt2 = 1, = 1.5 and /X4 = 2. 
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Fig. 3. — Distribution of the radiation power emitted at all frequencies F(fi) = J F(v,^)dv (lower panel) 
and at various frequencies F(fj,/x) (upper panel) over the dimensionless angle \x — k z cT/uj, where v\ — 0.5, 
v-i = 1, iaj = 1.5, V4 = 2 and v§ = 4. 



